Dirac Equation on a Curved 2 + 1 Dimensional Hypersurface 



Mehmet Ali OlpakB 

Department of Physics, Middle East Technical University, 06531, Balgat, Ankara, Turkey 

(Dated: December 26, 2011) 

Interest on 2 + 1 dimensional electron systems has increased considerably after the realization 
of novel properties of graphene sheets, in which the behaviour of electrons is effectively described 
by relativistic equations. Having this fact in mind, the following problem is studied in this work: 
when a spin 1/2 particle is constrained to move on a curved surface, is it possible to describe this 
particle without giving reference to the dimensions external to the surface? As a special case of this, 
a relativistic spin 1 /2 particle which is constrained to move on a 2 + 1 dimensional hypersurface of 
the 3+1 dimensional Minkowskian spacetime is considered, and an effective Dirac equation for this 
particle is derived using the so-called thin layer method. Some of the results are compared with 
those obtained in a previous work by M. Burgess and B. Jensen. 

PACS numbers: 



I. INTRODUCTION 

Graphene sheets, which are effectively 2 + 1 dimen- 
sional objects, have some interesting properties that are 
relevant to both experimental studies and theoretical in- 
terest. The electrons in graphene exhibit a linear energy 
momentum dispersion as if they were massless Dirac par- 
ticles [l[. Some experimental studies on graphene has 
been awarded the 2010 Nobel Prize in Physics @, and 
this further encourages both experimentalists and theo- 
rists in studying various aspects of effectively 2+1 dimen- 
sional electron systems. One can also find field theoretic 
approaches to graphene in the literature @ . In this work, 
a 2 + 1 dimensional model system with an electrostatic 
backgroud is considered and the dimensionality of the 
system is undertood as constraining a 3 + 1 dimensional 
system onto a 2 + 1 dimensional curved hypersurface by 
making use of the electrostatic background, as in 

In classical mechanics, there are certain types of prob- 
lems in which geometrical constraints (which one may as- 
sume to be holonomic for present purposes) are imposed 
on physical systems. When dealing with such problems, 
the nature of interactions that generate those constraints 
are not generally taken into account directly, since the 
problems can be handled by other means (generally us- 
ing the method of Lagrange multipliers). However, one 
does not deal with such cases in conventional problem 
sets of quantum mechanics, if, of course, one does not 
pay special attention to the issue. Beginning with Dirac, 
people tried to transpose that sort of problems to the 
framework of quantum mechanics (for example @, 
Q). Dirac developed a quantization procedure which 
involved handling the constraints within the commuta- 
tion relations @. Other relevant works involved 
making explicit use of the geometrical relations arising 
from the existence of the constraints @, ||, 0, named 
as the "thin layer" method by some authors [10| . How- 
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ever, the results obtained using those two approaches do 
not necessarily match with each other; and this is what 
one observes in the literature (see @, [Til] ) - There may 
be several possible solutions to this contradiction, which 
should eventually be verified experimentally. In any case, 
both approaches should be applied to a wider range of 
problems in order to be able to develop new ideas from 
the consequences. 

In this work, the above mentioned geometrical ap- 
proach will be applied to Dirac equation in 3 + 1 di- 
mensional Minkowskian spacetime M 4 to reduce it to a 
2+1 equation, which will be viewed as if written on a 
curved hyper surface of M . The problem has been previ- 
ously considered by M. Burgess and B. Jensen whose 
results are partly discussed in this work. The reason 
for choosing this approach is that it relies on the equa- 
tion itself, without searching for a more general quanti- 
zation procedure, and is thus only a limit of the 3 + 1 
dimensional equation. This makes it slightly more reli- 
able in the sense that if it succeeds in the task given, 
then an appropriate (effective) Lagrangian formulation 
and a quantization procedure may be developed begin- 
ning from the resulting equation. However, the nature 
of possible interactions/mechanisms that may generate 
geometrical constraints are not studied in detail, and a 
simple elctrostatic background is assumed as in [4| . The 
reader will notice that results presented in this work do 
not exactly coincide with the counterparts given in Q ■ 



II. GEOMETRY 

Although the problem deserves a general treatment, 
we will choose a simple and specific geometry, which was 
also used by da Costa in [(| and Ferrari and Cuoghi in 
0] • We refer the reader to Mitchell's treatment [l2j which 
covers possibly the widest range of different geometries. 

The above mentioned treatments [f| 0] consider non- 
relativistic particles confined to move on a 2 dimensional 
surface in 3 dimensional Euclidean space, and the ap- 
proach can be used to handle N — 1 dimensional hyper- 
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surfaces of the N dimensional Euclidean space directly. 
The geometry is expressed by the metric tensor trans- 
formed to a new basis via a general coordinate trans- 
formation. The curvilinear coordinates involved are ex- 
pected to have no special properties, except that one of 
the coordinates is chosen to be orthogonal to the surface 
on which the particle will be constrained. Then, the rel- 
evant Schrddinger equation is written in terms of these 
new coordinates, and the normal coordinate is squeezed 
so that one obtains an effective equation which involves 
only the surface coordinates. More explicitly 
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Here, r is the position vector of a generic point P on 
the surface, R is the position vector of a point Q just 
above the surface, N is the unit normal of the surface 
at P, are the curvilinear coordinates and the inner 
product is taken using the Euclidean metric 5 a i,, a,b = 
1,2,3. Notice that this expression is indeed exact with 
the coordinates chosen, but one may well write a Taylor 
expansion of position vectors [l3|. However, this task is a 
little bit more complicated, and for the present purposes, 
the above expression is perfectly useful. 

In these coordinates, one obtains the following normal- 
ization integral for the wave function [f| 0, HH : 



<FxV*(X)V(X) = j d 3 qVG^*(q)^(q) 

d 3 q^ (1 + q 3 Tr(a) + {q 3 f det(a)) tt*(ff)*(<?) = 1, 

(2) 



where X a are the Cartesian coordinates, G is the deter- 
minant of G/j, V) g is the determinant of the metric tensor 
dr dr 

oaa = — — ■ — — r induced on the surface, and a is the Wein- 
y 3 dq l df 

garten matrix of the surface 0] which is related to the 
extrinsic curvature of it and denned via @, 0, Ell : 



<9N 



dr 



(3) 



Here, and in the following parts, Einstein summation 
convention is used. 

The redefinition of the wave function within the nor- 
malization integral is crucial, because the factor coming 
in front of the new wave function carries the signature of 
the "external world" when substituted into the equation. 
With this procedure, after taking the limit q 3 — > 0, one 
may calculate any observable without giving any refer- 
ence to the external world [13| . 

In the treatment for Dirac equation, this geometry will 
be assumed to correspond to the spatial part of the space- 
time, and will be taken as constant in time, therefore 



causing no coupling between temporal and spatial parts. 
The only difference will be that Dirac equation is first 
order in all derivatives, while Schrddinger equation is sec- 
ond order in spatial derivatives, so we will use the above 
objects up to first order in g 3 , which also means there 
will survive no exact expressions for the metric tensor 
and related quantities. Such an expansion will do the 
job. 



III. DIRAC EQUATION 

We consider an electron (or a spin half particle) in flat 
spacetime, but the use of general curvilinear coordinates 
requires writing the Dirac equation as if a curved space- 
time is involved. This well known equation is written in 
curved spacetime in the following form [14], [13]: 

(i~f a E a - m)4> = 0, (4) 

where 7° are the well known Dirac matrices satisfying: 

[ 1 a , 1 b ] + = 2r l ab (5) 

plus denoting the anti-commutator and E a M are known 
as the inverse vierbeins satisfying: 

= E a »E b v r, ab , M , v,a,b = 0...3, (6) 

and is the appropriate covariant derivative which will 
explicitly be given below. Although the (inverse) vier- 
beins can be thought as coordinate transformation coef- 
ficients, this is not necessary in general. Indeed, these ob- 
jects appear in the equation due to the requirement that 
Dirac equation should be written using an orthonormal 
basis of vectors or one-forms in order to be able to use 
the flat spacetime gamma matrices [l4| . In addition, they 
are defined up to a local Lorentz transformation which 
leaves the relevant action invariant [l4j . 

In order to apply the thin layer method to this equa- 
tion, one needs to expand the included objects in powers 
of q 3 . This was performed before by Burgess and Jensen 
in [4|, but for surfaces having zero intrinsic curvature. 
When the surface has zero intrinsic curvature, one has 
the chance to use Cartesian coordinates in the vicinity 
of the surface [3]. However, this does not imply that the 
connection coefficients do not contribute to the geometric 
terms, contradicting with the statement of the authors. 
Indeed, making use of the algebra satisfied by the Dirac 
matrices, one observes that: 



due to 



1 C E C ^ aV [ 7 a , l b ] = 2 7 a (V^ Q ") (7) 



l a [l\ 71 = [7 C , l a ]l b ~ W + ^ hc l a (8) 



regardless of dimension of the spacetime. 

Now, let us write the equation explicitly. First, it 
should be noted that, the vierbeins and inverses satisfy 
the following (l4| : 



(9) 
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and the covariant derivative is given as [l4 |. [13j: 
Dp = dp +lo^, 
Up. = ^ QbA1 [7 a ,7 & ]; 

u a b p = -E b ^V M e a „, (10) 

where V is the covariant derivative in the so called co- 
ordinate basis which is not necessarily orthonormal, 
in other words: 

v l sr = d v y v + T v llX v x . (11) 

Now, we introduce the following set of vierbeins 
[lH (somehow different than in the reference): 

e a i = d lX a + q 3 8 t n a = d t x a + q 3 Oi k d k x a , 

e° o = 1, e° 3 = N a , a = 1, 2, 3, others = 0, (12) 

whose inverses are given by (to first order in q 3 ): 

E a % — E a 1 — q 3 a k l E a k , E = 1, E a 3 = N a , 
a = 1,2,3, others = 0, (13) 

where the over tilde implies that the quantity is evalu- 
ated on the surface and x a are the flat coordinates of 
the generic point P lying on the surface. This notation 
implies: 



~e a a = (dpX a ) . 



p — ywp.L ;u - (14) 

We also redefine the spinor in the following way [3] , [H( : 



X = ^^/l + q 3 Tr(a). 



(15) 



where the geometrical relations give Tr(a) = T\ 3 , i = 
1,2. Using these objects, one obtains the following equa- 
tion near q 3 — 0: 



of the spinor and the other from the connection coeffi- 
cients. This term is the analogue of the geometric poten- 
tial which appears in the Schrodinger equation in [a, Lfl] ; 
however, while that geometric potential was a scalar func- 
tion, this seems like a vector potential, but one which 
does not have a temporal component. This seems weird, 
of course; in the case of a specified interaction, the mean- 
ing of this term may become more clear. 

Burgess and Jensen assume that their resulting equa- 
tion is separable, and they consider two equations in 
one involving only the normal coordinate, and the other 
involving surface parameters. They conclude that the so- 
lutions behave like a shifted Gaussian along the normal 
direction. Though we proceeded in a somewhat different 
way, we may give an argument which is also in accor- 
dance with that of Burgess and Jensen; that is, we may 
assert that the solution to the 3 + 1 dimensional equation 
in the same coordinates should have an extremum, more 
precisely a maximum at q 3 — 0, if it is really possible to 
constrain the particle to that surface. One may also take 
this argument in the following way: X (t, r)\(t, r) is the 
probability density for the particle to be found at time t 
at the point r (remember that the spinor has been rede- 
fined so that the measure of the normalization integral 
involves only the determinant of the metric induced on 
the surface), and so, if the particle is constrained on a 
surface, then "it is most probably on the surface". So, 
the probability density should be maximum at q 3 = 0. 
These statements do not directly imply that x has a max- 
imum at q 3 = 0; however, it will be observed that a 
Gaussian distribution around the surface necessarily re- 
sults in that. To see this, let us introduce an electrostatic 
potential as a background, and derive the consequences: 



i h°do + l A E A i d i + - 1 A {V i E A *) 



1 



-y A N A Tr(a)+i A N A d 3 \ x 



(e 7 °^o(<Z 3 ) - m) X = 0, A = 1, 2,3; i = 1, 2, (18) 



- p A N A Tr(a) X + l A N A d^ - m* = 0, 
A = 1,2,3; i= 1,2, 



(16) 



where the over tildes imply the object is evaluated at 
q 3 = 0, and x c are the flat coordinates as introduced 
before. In order to obtain this equation, one makes use 
of (El) and: 



N a E a 1 = 0. 



(17) 



Here, the term — — r ) A N A Tr{a)x appears due to the 

existence of an external world, that is, it is the residue of 
the constrained dimension within the equation, there are 
two contributions to this term: one from the redefiniton 



where e is the charge of the particle. Following the argu- 
ments of [ij], one can identify this electrostatic potential 
with some A times q 3 , and expect that the behaviour of 
the spinor along the normal coordinate to be determined 
by the relevant derivative and the coupling term: 



(7W - i~f A N A d 3 ) X = K X . 



(19) 



Up to this point, every term but the normal dervative 
has been evaluated at q 3 = 0. However, there is an order 
q 3 term now. Here, one has to assume that eXq 3 also 
has an 0(1) contribution, which is consistent with the 
physics of the problem j4|. But there is another issue 
here: components of the normal vector are functions of 
the surface coordinates, so how can one get rid of the 
dependence on these coordinates in this equation? One 
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can simply expand the spinor in terms of the eigenvec- 
tors of the matrix Na and assume seperability of the 
equation for each term in the expansion. The eigenval- 
ues of this matrix are ±1, and assuming a Gaussian- like 
dependence for all componentes of the spinor imply: 



X = exp[-a(q ) + bq 



|X>(*> 

r=l 



<f)W), (20) 



where a > and V r (t, q l ) are the eigenvectors of i^ a Na- 
If the equation is seperable, K must be a constant ma- 
trix. It is a chance that one does not have to be that 
strict at all: K can be a function of time and the surface 
coordinates only. Putting this expression into (|19[) , one 
obtains: 

4 

exp[-a(q 3 ) 2 + bq 3 }J2[^l° 

r=l 

+ c r (~2aq 3 + b) - kr] K{t, q l )V r {q l ) = 0, (21) 

where c r and k r are the eigenvalues and corresponing 
seperation constants (or functions) respectively. Impos- 
ing that the probability distribution has a maximum 



at q a 



amounts to b = 0, which also amounts to 



k r = 0, Vr. This contradicts with the result in since 
in that reference the Gaussian distribution is shifted and 
the effect is interpreted as a contribution to the effective 
mass of the states after constraining the system. The 

solution is then a pure Gaussian with a = . The 

remaining (tangential) equation is then: 



E 



( 7 a + 7 A ^ % + \i a {sSa *)) 

(h r {t,q l )V r {q 1 )) =0. (22) 



— Tr(a) — m 
4 



The geometric term can still be interpreted as an ef- 
fective contribution to mass [H , and the equation can be 
written as a system of coupled 2x2 equations for two 
2-spinors. One expects that these equations should be 
decoupled for two linear combinations of those 2-spinors. 
Consider, for example, the following 4 representation for 
the gamma matrices: 

i r n .-_A1 

(23) 

where all entries are 2x2 matrices and <j a , A = 1,2,3 
are the well known Pauli matrices satisfying: 

[* A ,a B } + = 26 AB . (24) 

In this representation, after multiplying (|22p with 



"l " 







ia A 


-1 


,7 A = 








ij Na, the equations appear as the following: 

i( - c r a A N A d + a A E A l d z + ^a A (V z E A l ))v r {t, q l ) 

+ (jIY(a) +m^u r (t,q l ) =0, 
i(c r a A N A d + a A E A l d, + ^a A (ViE A % ))u r {t,q l ) 

+ (jTr(a) +rnjv r (t,q l ) = 0, (25) 



where h r (t,q l )V r (q l 



Though the equa- 



u r {t,q l ) 
[vr{t,q l )\ 

tions are still coupled, this representation in 3 + 1 dimen- 
sions naturally suggest the follwing representation for the 
(coordinate dependent) gamma matrices in 2 + 1 dimen- 
sions: 



±a A N A , T l 



±ia A E A 



A = 1,2,3, 



which clearly satisfy 

[r a ,r' 3 ] + = 2 ff ^, 



1,2, 
(26) 



(27) 



and g a P is clearly the induced metric on the 2 + 1 dimen- 
sional hypersurface. 

Here, there is one interesting fact: the charge of the 
particle does not appear in the reduced equations! This 
is interesting, because it suggests the idea that the elec- 
trostatic interaction in 3 + 1 dimensions seems to mimic 
gravity in 2+1 dimensions, since the considered geometry 
is curved and the extrinsic contributions can be absorbed 
into the mass term. This fact then brings a new ques- 
tion: how can one incorporate gravitation effectively in 
this picture? For sure, the answer will involve the com- 
plete set of conditions on the nature of the solutions to 
the tangential equation. 

However, in the case of a magnetic interaction, one may 
not be able to reduce the equation to 2 + 1 dimensions. 
This is suggested even by non-relativistic spin-magnetic 
field interaction: 



H, 



2mc 



(V x A • S) 



(28) 



where A is the electromagnetic vector potential. The 
cross product counts the number of spacetime dimen- 
sions, and a possible reduction scheme will most proba- 
bly be gauge dependent [l3[. Beyond speculations, one 
has to check what happens in this case. 



IV. CONCLUSION 

In this work, the treatment of Burgess and Jensen has 
been generalized to surfaces having nonzero intrinsic cur- 
vature, and the arguments on the origin of the geometric 
term and the behaviour of the solutions along the normal 
direction have been corrected. However, there still are a 
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number of unanswered questions for this case. First of all, 
solutions to the reduced equations should be derived and 
interpreted in a future work, and limitations on possible 
experimental realizations should be discussed. Secondly, 
one has to examine whether there are any geometrical 
constraints on the relevant surface, as done by Burgess 
and Jensen in their treatment. The last, and may be the 
most interesting unanswered question is, as indicated in 
the previous section, is the one concerning similarity with 
2 + 1 gravity. 

One can think of more interesting but complicated sce- 
narios in terms of the interaction constraining the system 
into 2 + 1 dimensions. However, the treatments seem to 
be limited in the case of single particle equations. A 
more general, may be a field theoretic treatment can be 



studied in order to develop deeper, and experimentally 
verifiable ideas. 
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